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Status of the Hadronic r Decay Determination of \Vus\ 
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We update the hadronic r determination of \Vua\, showing that current strange branching fractions produce 
results 2 — 3a lower than 3-family unitarity expectations. Issues related to the size of theoretical uncertainties 
and results from an alternate, mixed r-electroproduction sum rule determination are also considered. 



1. Introduction and Background 

The determination of \Vus\ from hadronic r 
decay data rests on the finite energy sum rule 
(FESR) relation, 

valid for any analytic w{s) and kinematic- 
singularity-free correlator, H, having spectral 
function, p{s). To obtain |T4s|, Eq- HI) is applied 
to the flavor-breaking (FB) correlator difference 



w{s) p{s) ds 
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where 



n 



are the spin J = 0, 1 components of the 
flavor ij, vector (V) or axial vector (A) current 
two-point functions, and (0 + 1) denotes the sum 
of ■/ = and 1 components. The OPE is to be 
employed on the RHS for sufficiently large sq- 

The spectral functions, Py/j^.^j, are related 
to the differential distributions, dRv/A\ijlds, 
of the normalized flavor ij V or A current 
induced decay widths, Rv/A-.ij = r[T^ 
UrhadTonsv/A-ij (7)]/r[r" UrerV^i'^)], by 
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matrix element, and, with Sew a short-distance 
electroweak correction [2], = 12tt'^ Sew /'m'i- 
Use of the J = -I- 1, FB difference An.^, 
rather than the analogous difference involving 
the linear combination of J = 0, 4- 1 spectral 
functions appearing in Eq. ([2]), is a consequence 
of the extremely bad behavior of the integrated 
J = (longitudinal) D = 2 OPE series [3]. 
Fortunately, apart from the accurately known tt 
and K pole terms, contributions to Py^^^.^ are 
oc [{rrii + TOj)^], making ud continuum contribu- 
tions negligible. Once the small continuum us 
J = contributions are determined phenomeno- 
logically using dispersive [4 and sum rule [5] anal- 
yses of the strange scalar and pseudoscalar chan- 
nels, respectively, the J — Q contributions can be 
subtracted, bin-by-bin, from dRy+A-ij /ds, allow- 
ing one to construct the re-weighted J = + I 
spectral integrals, Ry+A-iji^o), defined by 



Rv+A■,^Ji■^0) _ r 



dsw{s) p^y^2ij(^) ' (3) 



and, from these, the FB combinations, 

Rv+A;udi^o) Ry+A-usi'^o) 



SR^+A^o) = 
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dsw{s)Apr{s) . (4) 







With I Kid I, and any parameters in 
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from other sources, Eq. ([T]) yields [5] 



lK.sl = 



\ 



R 



V+A;u8 



(so) 



;ud(«0) 



(5) 



S-riW.OPE, \ 

on,^j_A [So) 



^V+A-ud,usi^o) 



v+A yoj is typically << 
(usually at the few-to-several-% level) for sq > 
2 GeV^, making a high precision |Kts| determina- 
tion possible with only modest OPE precision [6J . 

It turns out (see also below) that the conver- 
gence of the integrated J = 0+1, D — 2 OPE se- 
ries may also be somewhat problematic. As a re- 
sult, it is also of interest to consider FESRs based 
on the alternate FB correlator difference, 
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where Hem is the scalar part of the electromag- 
netic (EM) current two-point function. All a/ 
shares with AIIt- the vanishing of 13 = con- 
tributions to all orders but, by construction, has 
strongly suppressed D = 2 contributions [7]. 
D — A contributions turn out also strongly sup- 
pressed compared to those of AIIt-. This sup- 
pression does not, however, persist beyond D — 
4 [7]. The EM spectral function, pem{s), re- 
quired on the LHS of the AIIm FESR, is given by 
Pem{s) = j"^o(^) ^ with (7o(s) the bare inclusive 
hadronic electroproduction cross-section. The 
AIIm FESR yields a solution for \ Vus \ of the form 
Eq. (ini), with the RHS denominator replaced by 



where RemM 



^-^2,^(^0) — 



dsw{s)pEM{s) and 
dsw{s)AnM{s) . (7) 



'\s\=sn 

2. Spectral and OPE Input 

2.1. Spectral Input 

We compute Ry+A-udi^o) and Ry+A-usiso) us- 
ing the publicly available ALEPH us [S] and 
ud [9] spectral data and covariances. Separate 
ud V and A analogues, Ry/A-ud^^o), required for 
the mixed r-electroproduction FESRs implement 
the improved KKtt V/A ud separation [TIT pro- 
vided by CVC and the BaBar determination of 
7 = 1 KKtt electroproduction cross-sections [TT] . 



A small global rescaling of the continuum ud 
V-|~A distribution accounts for recent changes in 
Sew, Rv+a-us and Be- We employ |Kd| = 
0.97425(23) [Tg and current values 13] for B^, 
Rv+A;us and Rv+A;ud- Since BaBar and Belle 
have not yet completed their remeasurements of 
dRv+A:us/ds, we work with an interim partial 
update obtained by rescaling the 1999 ALEPH 
distribution [8] mode-by-mode by the ratio of 
new to old world averages for the branching frac- 
tions [14- The new world averages, based on the 
results of Refs. [15116117118119120121122) . are given 
in Table [1] [13j . The us V+A covariance matrix 
cannot yet be analogously updated, so the im- 
proved precision on the us branching fractions 
translates into an improved us spectral integral 



error only for w = w^^^rp = i(j(oo) and sq 

Details of the treatment of the EM spectral 
data, required for the spectral integral side of the 
AHm FESR, are omitted here because of space 
constraints, but may be found in Ref. [7j. 



to: 



Table 1 

Current world averages, B[t^ - 
main strange hadronic modes X 
where required. Column 3 gives the references 
associated with post-PDG2006 changes. 
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with S factors 
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0.290(18) {S = 
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K-T] 
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K-uj 
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(KAtt)- (est'd) 


0.011(7) 






K*~r] 


0.014(1) 
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K-cj) 


0.0037(3) {S = 
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[16|19] 


TOTAL 


2.845(69) 
(2.870(68)) 
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Figure 1. \Vus\ versus sq from the AII^ FESRs 
for, from top to bottom, W20, Wio, wiq and W(oo)- 
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2.2. OPE input 

To keep OPE-breaking contributions from the 
vicinity of the timehke point on \s\ = sq suffi- 
ciently suppressed, we restrict our attention to 
w{s) having at least a double zero at s — sq, and 
to So > 2 GeV2 [ig. 

The leading, D = 2, OPE contribution to AII^ 
is known to 4 loops |24j : 



3 rUsiQ^) 
27r2 



2 27r2 Q2 
19.93a2 + 208.75a3 + 



7 

l + -a 



(8) 



with a = as{Q^)/n, and as{Q^) and ms{Q^) the 
running coupling and strange quark mass in the 
MS scheme. Since independent determinations of 
as imply a{m'^) ~ 0.1, convergence at the space- 
like point on \s\ = so is marginal at best. With 
such slow convergence, conventional prescriptions 
for assessing the D = 2 truncation uncertainty 
may lead to significant underestimates. 

To deal with the potential D ^ 2 convergence 
problem, one may either work with AUt and 
w{s) chosen to emphasize regions of the com- 
plex s = —Q^-plane away from the spacelike 



point, where |q;s(Q^)| is smaller and convergence 
improved pS], or switch to the alternate AIIm 
FESRs where D = 2 contributions are suppressed 
already at the correlator level [7j. In the latter 
case, the D = 2 contribution becomes 

, 2xiOPE 3 m^(g2) 



+ 4.384a2 



27r2 
UMa^ 4 



1 



(9) 



more than an order of magnitude smaller than in 
the An^ case. Since as(so) grows with decreas- 
ing So, making higher order terms relatively more 
important at lower scales, extracted \Vus\ results 
will display an unphysical sp-dependence if ne- 
glected, higher order D = 2 terms are, in fact, 
not negligible, so-stability studies thus provide 
a handle on the impact of the potentially slow 
integrated D = 2 convergence on \Vus\- 

D = 4: OPE contributions to An^(Q2) and 
AHm {Q'^) are determined by (m^ss) and {mii£), 
up to negligible O(to^) corrections. The rele- 
vant expressions, as well as those for the D = 6 
four-quark condensate contributions, are easily 
constructed from the results of Ref. [26], and 
given in Ref. 7J. If one works with weights 
w{s) = Z]m=o^™y™' ■^it^^ y = ^/^o, integrated 
D = 2k + 2 OPE terms scale as I/sq, allowing 
contributions of different D to be distinguished 
by their differing so-dependences. 

As £) = 2 OPE input, we employ ms{2 GeV) = 
96(10) MeV dZj and asiml) = 0.323(9), the 
latter obtained from an average, as{Mz) = 
0.1190(10), of various recent determinations (in- 
cluding lattice 28J and r [29] results, which are 
now in very good agreement) via the standard 
combination of 4-loop running and 3-loop match- 
ing at the flavor thresholds [5U] . 

At D = 4, we employ the GMOR relation for 
{niiU) and evaluate (m^ss) using the ChPT de- 
termination ofnis/mi [3r and (msSs) / {mi££) = 
1.2(3), the latter obtained by updating Ref. ^2] 
using the average of recent Uf = 2-1-1 lattice 
determinations of fsJIs as input [33] . 

D = Q contributions are estimated using 
the vacuum saturation approximation (VSA), 
rescaled by pvsA = 1(5), while D > Q contribu- 
tions are neglected. Since integrated _D > 6 OPE 
contributions scale as 1/s^ [N > 2), if > 4 
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Figure 2. \Vus\ versus sq from the An^- FESRs 
for, from top to bottom, W20, i&io, wio and W(oo)7 
with the spectral input modified by rescahng up 
by 3cr the branching fraction of the large, but not 
yet remeasured, K^t:~'i:^ mode. 



Figure 3. |Kis| versus so from the AE^ FESRs 
for, from top to bottom, W20, wiO: ^iq and i/;(oo), 
with the spectral input modified by rescaling up 
by 3(7 the branching fractions of all modes not yet 
remeasured by either BaBar or Belle. 
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contributions are, in fact, not small, and these 
input assumptions are unreliable, an unphysical 
so-dependence of \Vus\ will result, again making 
so-stability tests important. 

3. Results and discussion 

The results for | Vus \ obtained using the inputs 
specified above for the AIIt- FESRs based on the 
J = + 1 kinematic weight W(oo) [y] , and three 
weights, wio(y), W2o{y), wioiy), constructed in 
Ref. [25] specifically to improve the poor inte- 
grated J = 0+1, D = 2 convergence, are dis- 
played in Fig. [TJ The so-instability of the W(oo) 
results is much greater than the theoretical uncer- 
tainty ^ ±0.0005 often quoted for the sq = 
version of this analysis in the literature. The re- 
sults corresponding to wio, in contrast, display 
a very good window of SQ-stability. A positive 
feature of the AH,- analysis is the fact that the 
results for all four weights appear to be converg- 



ing towards the stable wia value as sq — > . The 
So — rn^ versions of the various analyses are 



IK, 



0.2180(32)(15) (wio) 

0.2188(29)(22) (w2o) 

0.2172(34)(11) (wio) 

0.2160(26)(8) (w;(oo)) 



(10) 



where the first error is experimental (dominated 
by us spectral errors) and the second the nom- 
inal theoretical error. The nominal theory error 
is obviously much smaller than the observed so- 
instability in the W(oo) case, and hence unrealis- 
tically small. Comparison to the results of earlier 
An^ FESR analyses [6,2134 35j shows the signif- 
icant impact of recent, improved us experimental 
results on the \Vus \ central values. The decreases 
represented by the remeasured us branching frac- 
tions, lead to \Vus\ results 2 — 3(t below the 3- 
family-unitarity expectation, 0.2255(1) [12] . 

It should be stressed that several important 
strange decay modes have yet to be remeasured 
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by either BaBar or Belle, and that the level 
of consistency of the Sq = results for dif- 
ferent weights could be significantly affected by 
such future remeasurements. As an illustration, 
we show, in Figure [2l the impact on \Vus\ as 
a function of sq of rescaling upward by 3(t the 
as-yet-unremeasured K'^tt^tt'^ branching fraction, 
and hence also the X'^tt^tt^ component of the 
us spectral distribution employed above. The 
issue of whether plausible shifts in the as-yet- 
unremeasured branching fractions are capable of 
restoring agreement with 3-family-unitarity ex- 
pectations is less clear. In fact, it would take si- 
multaneous 3(7 upward rescalings of all currently 
unremeasured us branching fractions to restore 
agreement. Such a rescaling, moreover, does not 
produce a convincing so-stability plateau for any 
of the weights considered, as shown in Figured 
The results for |Ks | for the AHm FESRs based 
on W(oo)(y), the weight wio{y) displaying the best 
So-stability for the AXI^r FESR, and the weights 

W2, W3 and W4, where wn = 1 - t^zjU + iv^' 
are displayed in Fig. 3] The weight wn produces 
a single surviving integrated D — 2N + 2 > 4 
OPE contribution suppressed by the coefficient 
1/{N — 1) and scaling as making it a useful 

choice in this case, where the slow integrated D = 
2 convergence found for the wn versions of the 
AIIt- FESRs is no longer relevant. 

If it was poor D — 2 convergence which was re- 
sponsible for the SQ-instability of the u'(oo) AII^ 
FESR results, one would expect to see a much 
improved stability plateau for the corresponding 
AIIm FESR, as is indeed found. The very good 
stability for the wn results also indicates that the 
integrated D — 2N + 2 contributions relevant to 
these cases become negligible in the upper part 
of the So window displayed in the Figure. Since, 
however, D > 6 contributions increase in going 
from AIIt- to AIIm [7], one would expect the in- 
stability for weights like liio, which do not sup- 
press these to the same extent as do the other 
weights considered, to be enhanced, as is indeed 
found to be the case. Even so, the Wio results 
converge well to the stable results for the other 
weights as so m^. 



Figure 4. \Vus \ as a function of sq from the mixed 
r-electroproduction FESRs for, from top to bot- 
tom at the left, wio, li'(oo), w^, W4 and W2- 
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Given the very good stability of the u'(oo) re- 
sults, it is possible to quote a final result based 
on the So = version of the u;(oo) FESR, which 
allows us to take advantage of the improvements 
in the us branching fraction errors. The result is 

IK.I = 0.2208(27)(28)(5)(2) (11) 

where the first three errors are due to the uncer- 
tainties on the us V-l-A, residual / = EM and 
residual ud V/A spectral integrals, respectively, 
and the fourth is due to the D = 2 and 4 OPE 
uncertainties (see Ref. [7] for further details). 

We conclude by stressing that, for both the 
AIIt and AIIm FESRs, improved errors on 
dRv+A-us/ds are crucial. This requires both 
remeasurements of as-yet-unremeasured strange 
mode branching fractions, pursuit of higher mul- 
tiplicity modes with branching fractions down to 
the fewxlO^^ level, and, in particular, a full in- 
vestigation of the K3tt and KAtt modes, which 
were not in fact measured, but rather estimated, 
in the earlier experimental analyses. 
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